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The flurry of activity in non equilibrium statistical phenomena covers many fields of theoretical 
and practical importance such as growth models, front propagation, crack propagation and many 
more. In this article we derive an exact inequality relating the response function, measuring the 
steady state response of the physical field of interest to an external probe, and the correlation 
function. With no further assumptions, that inequality is turned into an exponent inequality. For 
growth models of the general KPZ family at steady state, the response and dynamic exponents are 
equal. This leads to a more stringent inequality that relates the roughness and growth exponents. 
To demonstrate the usefulness of the inequality, we consider the NKPZ family of growth models 
and show that various approaches to that problem produce exponents, which violate the inequality. 
We also consider some experimental and numerical results that also seem to violate the inequality 
and suggest measurements that would clarify the origin of the difficulty. 
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The focus of interest in statistical physics has shifted in 
the last two decades from equilibrium phase transitions 
and later the dynamics of phase transitions 1] to the 
study of non-equilibrium systems such as various growth 
models 0, 0, 
propagation 



pi, front propagations crack 
I, etc. In spite of that shift the 
main objects of study remained of a similar nature, a 
small set of exponents which describe the steady state 
properties as well as the evolution of the system. Except 
for a number of one dimensional exactly soluble prob- 
lems [2, ^, IB, [H, 13, Mii^), the sets of exponents given 




in the literature for many systems belonging to the above 
categories, vary considerably from author to author and 
depend strongly on the method of derivation. This is 
very different from the situation in equilibrium phase 
transitions, where methods as different as high tempera- 
ture expansion, momentum space RG and real space RG 
yield very close exponents. Under such circumstances 
rigorous results that can put bounds on the exponents 
describing the system are obviously most valuable. In 
the following we present a quite powerful inequality for 
dynamical stochastic systems, which is an extension of 
the Schwartz-Soffer inequality derived for quenched ran- 
dom problems [l^ . The inequality is of a generic nature 
and relates the response at steady state of some mea- 
surable physical field to an external disturbance to the 
time dependent correlations of that physical field. This 
in itself is enough to check approximation schemes or ex- 
periments that supply both quantities. In the interesting 
cases, where the system may be described in terms of set 
of exponents, the predictions of the inequality become 



more dramatic by turning it into an exponent inequality. 
In order to assess the usefulness of the exponent inequal- 
ity, we apply it, as an example, to the results obtained 
by different methods for one specific model, the non- 
local Kardar-Parisi-Zhang (KPZ) equation This test 
demonstrates not only the obvious discrepancies among 
results obtained by different methods, as discussed above 
but also that all methods but one, namely the self consis- 
tent expansion (SCE) violate the inequality 
at some region of parameter space. This means that the 
inequality is strong enough to point out problems in var- 
ious approximation schemes, which may be the first step 
leading either to improving or discarding those schemes. 
Finally, we consider a number of experimental and nu- 
merical studies that yield exponents in violation of the 
inequality and suggest measure- 
ments that would clarify the origin of the problem. 

Many interesting dynamical physical systems may be 
described in terms of some physical field, (j){r, t) driven by 
a "noise" field, ri{r,t). The fist of systems, described by 
generic Langevin field equations, includes models of crit- 
ical dynamics [H, growth models of the KPZ family Q 
and its many variants 0, S Q , noise driven Navier Stokes 



24| etc.. Strictly speaking, the physical field given as a 



function of time depends not only on the noise field at 
earlier times but also on initial conditions. The depen- 
dence on initial conditions decays, however, in time, and 
we are left with an implicit relation between the Fourier 
transform of the field and the Fourier transform of the 
noise, 
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(/)(q,cj) = 0{q,w; 77(1,(7)} , (1) 
where 77(1, a) is a Gaussian random field with {ri{l, a)) = 



2 



0, and 

(77(1, a)77(m, ^)) = 2Do{l, a)6{l + m)S{<j + . (2) 

We use the traditional convention and take the physical 
units of the noise to be equal to the physical units of the 
physical field over time. The loss of initial condition 
information is characteristic of the steady state and thus 
the various averages we will consider in the following, 
which are averages over the noise are also steady state 
averages, traditionally studied in the literature 

We are interested in the response fimction G{q,Lu), to 
be defined by 

(||^)-G(,,.)%-p)^(.-a), (3) 

and in the correlation function ^{q,Lu), which is defined 
by 

(0(q, t^)</)(-p, -a)) = $((?, uj)S{q - p)S{uj - a) . (4) 

Because of the Gaussian character of the noise, using 
integration by parts the response function can be written 
also as 

G(g,w)(5(q-p)(5(w-o-) = (0(q, w)7y(-p, -a)) /2Do{q,uj) . 

(5) 

Note, that if we define the response function by the 
right hand side of equation ([5]) (which still involves a 
non-trivial correlation function) the rest of the deriva- 
tion follows, even if the distribution of the noise is not 
Gaussian. 

The average (x(q, 'i^)'0("Pi ^cr)} can be viewed as a 
scalar product of x(q, and ^lJ{p,a), since it has all 
the properties required of a scalar product. Using the 
Schwartz inequality we find 

\G{q,cj)\S{<i-p)5{cj~a) (6) 
< ^^q,u;)2Do{q,u;)6{ci - p)Siu; - a)/2Doiq, iu) . 

Integrating over p and a and squaring both sides leads 
to 



2G(q, uj)Gi-q, -u)Do{q, u) < . 



(7) 



The above is a general exact inequality, relating the 
response function, as defined by the right hand side of 
equation ^ and the correlation function. To turn that 
into an exponent inequality, let the equal time correla- 
tion. 



A(g) 



(8) 



be proportional to q^^ for small q and let the response 
exponent z characterize the small q behavior of the re- 
sponse function 



The characteristic frequency, Ld{q), associated with the 
decay in time of the correlation is given in terms of the 
dynamic exponent z as [1] 



a>-i((Z)=^<&(q,0)/A(g) 



(10) 



and is proportional to q^ . (The use of a small q power 
law behavior of A, G{q, 0) and lo is not a restriction at 
all as any dependence that is not a power law at small q 
can be described as a power law with a power which is 
either zero or ± infinity.) 

We will concentrate in the following on bare spectral 
functions Do{q,u!) (the noise correlator in Eq. ([2])) that 
for small q and lo have the form 

Do{q, w) = Bq^^'^with a >0 for small q and cj . (11) 

The above form is rich enough to make our point. Still, 
the discussion of more general cases is straightforward. 

The required exponent inequality is obtained now by 
setting = in Eq. 



2z + 2<7 <T + z. 



(12) 



The response exponent z is usually related to other expo- 
nents characterizing the system. For example in Model 
A of Hohenberg and Halperin Q z = F. The latter is an 
exact relation, which leads to the correct but not very 
exciting inequality F < z. In growth models of the KPZ 
family, we have, instead, the relation, z = z, It 
does not seem, however, very likely that the two options 
mentioned above, z = T and z = z exhaust all the in- 
teresting possibilities. Independent determination of z in 
other systems might be thus of interest. 

Note that while inequality is rigorous, we will 
use from now on additional scaling relations, relating the 
exponents, which are most probably correct and widely 
used but have never been rigorously proven. Using the 
relation z = z, the general inequality above reduces to 
the simpler 



z<T -2(7 <T. 



(13) 



G{q,0)^q- 



(9) 



It is well known that for members of the KPZ family, 
there is a lower critical dimension above which appears a 
weak coupling solution on top of the strong coupling solu- 
tion, which is the only possibility at lower dimensions . 
The above inequality is correct for any set of exponents, 
whether strong or weak coupling and in fact for weak 
coupling exponents, we expect the leftmost inequality to 
hold as an equality. 

To demonstrate the usefulness of the exponent inequal- 
ity, we use it to test various theoretical approaches to the 
study of the KPZ family. We concentrate on the non- 
local KPZ, where we expect that there is a good chance 
that inherent problems in various theoretical approaches 
will show up clearly in violations of the inequality. We 
find that that is indeed the case and all the methods 
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17, 18, hJI violate the in- 



tested apart from the SCE 
equahty. 

The Non-local KPZ (NKPZ) equation has been intro- 
duced in Ref. [3| to account for the non-local hydrody- 
namic interactions in a deposition of colloidal particles 
in a fluid. It was later generalized to spatial correlated 
noise (cr ^ 0) in Ref. [5|. The NKPZ equation for the 
height function, h{r,t) is given by 



dthir,t) 



i^V^h+ J dr'g (r') V/i (r + r', t) ■ V/i (i 



where the kernel g {f) has a short range part, Ao(5'^ (r), 
and a long-range part ^ Xpr^^'^ . In Fourier space, g{q) — 
Ao -l- XpQ^'' ■ In order to avoid an extremely, complicated 
phase diagram, we discuss here only the case Aq = 0. 
The noise has zero mean, but is allowed to have spatially 
long range correlations, given by equation pip . (The 
roughness exponent, a, we are interested in, is related 
in KPZ to the exponent F by 2q! = F — d, where d is 
the dimensionality of the system.) The strong coupling 
solution found by the Dynamical Renormalization Group 
(DRG) is 0] 

^ , {d-2-2p){d-2~3p) 
"^^^ = ' + (3 + 2-P)d-6-9p ■ 

(Because of the extra scaling relation between a and z 
0, Q , namely a-\-z = 2 — p, there is only one independent 
exponent in NKPZ.) 

The above result violates the inequality over a 
whole range of parameters defined by T^jiq — z^j^q — 
2(7 < or explicitly given for the weaker inequality, with 
(7 = 0, by 

^ < P< Po{d) for d<dQ 



Po {d) < p< ^ for d> do 



with do ~ 3.395 and po {d) = ^ + 



(16) 



"din 2 



In 2 ' 

The shaded re- 



where W{x) is the Lambert function, 
gion in Fig. [1] is the region in the (d, p) plane where the 
weaker inequality (i.e. Eq. with cr = 0) is violated. 

In Fig. [2] we give the solution for the dynamic exponent 
z for d ~ 1, 2, 3 by an improved mode-coupling approach 
originally due to Colaiori and Moore for the local case 
[27i] . which was later applied by Hu and Tang 28 1 to 
the non-local case. The shaded area bounded by the line 
z = (d -f 4 — 2p) /3 marks the region where the inequality 
is violated. 

Note, that in both of the above examples the boundary 
of the violation is a independent. Namely, we have used 
the weaker inequality with a — 0, T > z. The use of the 
stronger inequality (Eq. (|13p with positive cr) puts more 
stringent bounds and therefore, the actual regions where 
violations occur are larger, depending on the value of a. 

Next, we check for violations of the inequality in results 
obtained by Tang and Ma J2d\ generalizing the Flory- 
type Scaling Approach (SA) of Hentschel and Family ^3Q] 



(14) 
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FIG. 1: Violation of the response-correlation inequality (|13p 
with cr = by the DRG method derived in Refs. 0, [El occurs 
in the shaded area enclosed by the curves p — (d — 2)/2 and 
p = po{d) in the phase diagram. 



to the non-local case. The strong-coupling dynamic ex- 
ponent obtained by using this method is 



ZSA 



(2-p)(d + 2-2cr) 
d + 3 - 2(T 



(17) 



It turns out that this solution violates the inequality in 
a whole region of parameter space defined by 



1 



d d 
2<"<^ 



1 



for p > — 1 . 



P d 

< cr < - for p < -1 . 

2 2 ' 



(18) 



(19) 



These results are presented graphically in Fig. [3] below. 

Last we test the results for NKPZ of the Schwartz- 
Edwards Self-Consistent Expansion (SCE) [13, El. This 



method predicts a whole zoo of possible phases summa- 
rized in Table I in Ref. 19] . It turns out that these results 
are not only consistent with the exact one-dimensional 
result fill but the exponents in all the phases are con- 
sistent with our stronger inequality over the whole pa- 
rameter space and all dimensionalities. In Fig. 3] below 
we give the SCE solution for the dynamic exponent z for 
d — 1,2,3 for uncorrelated noise (cr — 0)(but as stated 
above the inequality holds also for cr > 0). 

Our inequality, is shown thus to be a useful tool in 
detecting shortcomings of various analytical methods ap- 
plied to the NKPZ family and consequently suggest that 
results obtained by those methods should be suspected 
even in cases where the inequality is not violated, at least 
until the origin of violation of the inequality is under- 
stood. 

We have reviewed the literature and found some nu- 
merical and experimental results that violate the inequal- 
ity (fT3j) . We have found violations for crack surfaces in 
wood 3 and in mortar [l3| , where the dynamic ex- 
ponent z and the roughness exponent a are reported. In 
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FIG. 2: The Dynamic exponent z as a function of the non- 
locality parameter p] using the mode coupling result (blue 
line - reproduced from Ref. H^). The part of the blue line 
within the shaded region is excluded by the inequality (|13p 
with a = 0. 



addition, we have found violations in a number of sur- 
face growth systems, including interface roughness in Fe- 
Cr superlattices [20| . growth front roughening in sihcon 
nitride films by plasma enhanced chemical vapor deposi- 
tion 21 1, interface roughening in wrinkly metal (2^ and 
effects of shadowing on growth processes [2^, where all 



three exponents a, (3 and z are measured independently. 
The simplest explanation for this violation could be just 
due to one or more of the numerous artifacts in the meth- 
ods used to extract the exponents. Possible artifacts are 
discussed thoroughly in Refs. 31, 32, 3^. An alternate 



explanation for the problematic exponents could be that 
the time involved in the measurement is too short for the 
system to reach steady state while the roughness expo- 
nent a, should be measured in steady state. 

There are, however, two more alternatives. Exper- 
iments and numerical simulation often produce non- 
steady-state quantities. For example, 




The nonlocality parameter p 



FIG. 3: .Violation of the response-correlation inequality (|13p 
by the SA method derived in Ref. 29] occurs in the shaded 
area enclosed by the curves p = —1 and (j = (d+l + p)/2in 
the phase diagram. 



with h{v, 0) — 0. This tends to a steady state quantity 
only as time tends to infinity. Any q dependent charac- 
teristic time obtained from Eq. ([20]) . can be related to 
the steady state exponent z only if a single characteris- 
tic time exists for each length scale. In the presence of 
more than one characteristic time-scale for each length 
scale, the relation between steady-state and non-steady- 
state dynamic exponents is not obvious and our inequal- 
ity relating z to a cannot be invoked. Moreover, even for 
steady-state we have to recall that our inequality fT^]) 
is based on the relation z ~ z. For the KPZ family we 
understand the origin of that relation but for other sys- 
tems this might not be necessarily true (as for example in 
critical dynamics (H). Therefore only the more general 
inequality (dH), involving z and z, should be used. To 
test the exponents obtained by one method or another, 
one would need to obtain not only the time dependent 
steady state correlation function, 



(21) 



(22) 



$(r,0 = ([Mr,t)-MO,0)]^ 

but also the corresponding response function, 

5h {r,t)' 



G(r,t) 



(577(0,0) 



W{r,t) = ^[h{r,t) ~ h{0,t)Y 



(20) 



The fact that the scaling relation z = a/P is violated in 
the surface growth systems mentioned above might per- 
haps indicate the existence of more than a single char- 
acteristic time scale, but we cannot be confident of that. 
In fact it might also be an indication that steady state 
has not been reached. In order to clarify this issue we 
suggest that measurements such as those described in 
Refs. El, mill, would be complemented by steady 
state measurements of time dependent correlation func- 
tions and response functions given by equations (j2ip and 
([22|) . This does not require much additional effort, since 
in order to obtain a, the system has anyhow to be studied 
for times long enough to reach steady state. The inde- 
pendent extraction of the exponents z and z will enable 
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FIG. 4: (color online) The dynamic exponent 2; as a function 
of the non-locality parameter p for uncorrelated noise {a — 0) 
in d = 1,2,3 dimensions resulting from the Self-Consistent 
Expansion [Q. For d — 1 the solution saturates the inequal- 
ity for —1/2 <p<— 1/5. Note that for d = 2,3 two phases 
are possible in some range of the parameters as seen from 
the fact that z is a multivalued function of p (solid (blue) line 
weak coupling solution, and dashed (red) line strong coupling 
solution). Also note that for d = 2 and p = there is only 
one phase as the solid (blue) line does not exist for p = 0. 



to check (a) whether the non steady state z obtained for 
example from the correlation (|20p and the steady state z 
are equal and (b) if z equals z. This and the use of the 
more general inequality (jl2p will put the quality of the 
data to a much more stringent test. Violation of inequal- 
ity (|12p will then certainly indicate that the exponents 
obtained are inaccurate for one reason or another. 

To summarize, in this paper we showed how to gen- 
eralize the Schwartz-SofFer inequality derived originally 
for quenched random systems [l6| to dynamical stochas- 
tic systems. We show that the inequality, which involves 
the correlation and the response functions can be trans- 
lated into a simple inequality for the scaling exponents 
r, z and the noise correlation exponent a, in cases like 
the KPZ family, where z = z. Although being extremely 
simple, this inequality can be quite powerful when exam- 



ining analytical, numerical and experimental results. 

To demonstrate the utility of the inequality, we re- 
viewed analytical results for the non-local KPZ model, 
obtained by using four different methods: Dynami- 
cal Renormalization Group, Mode-Coupling, Scaling- 
Approach and the Self-Consistent Expansion. Interest- 
ingly, the first three methods yield results which contra- 
dict the inequality for a whole set of parameters, while 
the Self-Consistent Expansion is the only one which never 
violates it. This has an important implication on the 
choice of analytical tools when dealing with such stochas- 
tic models. 

Last, we reviewed a set of numerical and experimen- 
tal results coming from systems of surface growth and 
post-mortem analysis of crack surfaces. We suggest to 
measure the steady state response and correlation func- 
tions, which for surface growth systems, does not seem 
to require much additional effort. In case the more gen- 
eral inequality is violated it will become clear that the 
quality of the data and/or the methods of obtaining the 
exponents from it are not satisfactory. 

An interesting open question is whether there exist in- 
teresting cases where the response exponent z is not equal 
to r or to the dynamic exponent z. 

Hopefully, the simplicity yet strength of this result will 
also motivate researchers to explore the usefulness of rig- 
orous inequalities and derive improved ones alongside the 
more popular chase of approximated equalities. 
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